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STRESS- INTENSITY  FACTOR  EQUATIONS  FOR  CRACKS  IN  THREE-DIMENSIONAL 
FINITE  BODIES  SUBJECTED  TO  TENSION  AND  BENDING  LOADS 

1  2 
J.  C.  Neman,  Jr.  and  I.  S.  Raju 

SUMMARY 

Stress-Intensity  factor  equations  are  presented  for  an  embedded  ellipti¬ 
cal  crack,  a  semi-elliptical  surface  crack,  a  quarter-elliptical  corner  crack, 
a  seal-elliptical  surface  crack  along  the  bore  of  a  circular  hole,  and  a 
quarter-elliptical  corner  crack  at  the  edge  of  a  circular  hole  In  finite 
plates.  The  plates  were  subjected  to  either  reaote  tension  or  bending  loads. 
The  stress-intensity  factors  used  to  develop  these  equations  were  obtained 
from  previous  three-dlaenslonal  finite-element  analyses  of  these  crack  con¬ 
figurations.  The  equations  give  stress- Intensity  factors  as  a  function  of 
parametric  angle,  crack  depth,  crack  length,  plate  thickness,  and,  where 
applicable,  hole  radius.  The  ratio  of  crack  depth  to  plate  thickness  ranged 

from  0  to  1,  the  ratio  of  crack  depth  to  crack  length  ranged  from  0.2  to  2, 

and  the  ratio  of  hole  radius  to  plate  thickness  ranged  from  0.5  to  2.  The 

effects  of  plate  width  on  stress-intensity  variations  along  the  crack  front 

were  also  Included,  but  were  either  based  on  solutions  of  similar  configura¬ 
tions  or  based  on  engineering  estimates. 


fracture  strength,  accurate  stress-intensity  factor  solutions  are  needed  for 
these  crack  configurations.  But,  because  of  the  complexities  of  such  problems, 
exact  solutions  are  not  available.  Instead,  investigators  have  had  to  use 
approximate  analytical  methods,  experimental  methods,  or  engineering  estimates 
to  obtain  the  stress-intensity  factors. 

Very  few  exact  solutions  for  three-dimensional  cracked  bodies  are 
available  In  the  literature.  One  of  these,  an  elliptical  crack  in  an  infinite 
solid  subjected  to  uniform  tension,  was  derived  by  Irwin  [1]  using  an  exact 
stress  analysis  by  Green  and  Sneddon  [2].  Kassir  and  Slh  [3],  Shah  and 
Kobayashl  [4],  and  Vljayakumar  and  Atluri  (5]  have  obtained  closed-form 
solutions  for  an  elliptical  crack  in  an  infinite  solid  subjected  to  non- 
unlform  loadings. 

For  finite  bodies,  all  solutions  have  required  approximate  analytical 
methods.  For  a  semi-circular  surface  crack  in  a  semi-infinite  solid  and  a 
semi-elliptical  surface  crack  in  a  plate  of  finite  thickness.  Smith,  Emery, 
and  Kobayashl  [6],  and  Kobayashl  [7],  respectively,  used  the  alternating 
method  to  obtain  stress- intensity  factors  along  the  crack  front.  Ra  Ju  and 
Sewman  [8,9]  used  the  finite-element  method;  Hellot,  Labbens,  and  Pelllssier- 
Tanon  [10]  used  the  boundary-integral  equation  method;  and  Nlshloka  and 
Atluri  [ll|  used  the  finite-element  alternating  method  to  obtain  the  same 
information.  For  a  quarter-elliptic  corner  crack  in  a  plate,  Tracey  [12]  and 
Pickard  [13]  used  the  finite-element  method;  Kobayashl  and  Enetanya  [14]  used 
the  alternating  method.  Shah  [15]  estimated  the  stress-intensity  factors  for 
a  surface  crack  emanating  from  a  circular  hole.  For  a  single  corner  crack 
emanating  from  a  circular  hole  In  a  plate.  Smith  and  Kullgren  [16]  used  a 
f Inlte-element-alternatlng  method  to  obtain  the  stress-intensity  factors. 
Hechmer  and  Bloom  [17]  and  Ra  ju  and  Newman  [18]  used  the  finite-element  method 
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for  two  syaaetrlc  corner  creeks  emanating  froa  e  hole  In  e  plate.  Most  of 
these  results  were  for  Halted  ranges  of  paraaeters  and  were  presented  In  the 
fora  of  curves  or  tables.  For  ease  of  coaputation,  however,  results  expressed 
In  the  fora  of  equations  are  preferable. 

The  present  paper  presents  equations  for  the  stress-intensity  factors  for 
a  wide  variety  of  three-diaenslonal  crack  configurations  subjected  to  either 
unifora  reaote  tension  or  bending  loads  as  a  function  of  paraaetrlc  angle, 
crack  depth,  crack  length,  plate  thickness,  and  hole  radius  (where 
applicable);  for  exaaple,  see  Figure  1.  The  equations  for  uni fora  reaote 
tension  were  obtained  froa  Reference  19.  The  tension  equations,  however,  are 
repeated  here  for  coapleteneas  and  because  the  correction  factors  for  reaote 
bending  are  aodlf icationa  of  the  tension  equations.  The  crack  configurations 
considered,  shown  In  Figure  2,  Include:  an  eabedded  elliptical  crack,  a  seal- 
elltptlcal  surface  crack,  a  quarter-elliptical  corner  crack,  a  seal-elliptical 
surface  crack  at  a  circular  hole,  and  a  quarter-elliptical  corner  crack  at  a 
circular  hole  In  flnlte-thlckness  plates.  The  equations  were  based  on  stress- 
intensity  factors  obtained  froa  three-dlaens tonal  flnite-eleaent  analyses  [8, 
9,  18,  and  19]  that  cover  a  wide  range  of  configuration  paraaeters*  In  soae 
configurations,  the  range  of  the  equation  was  extended  by  using  stress- 
intensity  factor  solutions  for  a  through  crack  In  a  similar  configuration.  In 
these  equations,  the  ratio  of  crack  depth  to  plate  thickness  (a/t)  ranged  froa 
0  to  l,  the  ratio  of  crack  depth  to  crack  length  (a/c)  ranged  froa  0.2  to  2, 
and  the  ratio  of  hole  radius  to  plate  thickness  (r/t)  ranged  froa  0.5  to  2. 

The  effects  of  plate  width  (b)  on  stress-intensity  variations  along  the  crack 
front  were  also  Included,  but  were  either  based  on  solutions  of  similar  con¬ 
figurations  or  based  on  engineering  estlaates. 


NOMENCLATURE 


depth  of  crock 

width  or  half-width  of  cracked  plate  (eee  Fig.  2) 
half-length  of  crack 

boundary-correction  factor  for  corner  crack  in  a  plate  under  tension 

boundary-correction  factor  for  corner  crack  at  a  hole  in  a  plate  under 
tension 

boundary-correction  factor  for  embedded  crack  in  a  plate  under  tension 

boundary-correction  factor  on  stress  intensity  for  renote  tension 

boundary-correction  factor  for  surface  crack  in  a  plate  under  tension 

boundary-correction  factor  for  surface  crack  at  a  hole  in  a  plate 
under  tension 

finite-width  correction  factor 

angular  function  derived  froa  embedded  elliptical  crack  solution 

curve  fittLng  functions  defined  in  text 

bending  nultlpller  for  corner  crack  in  a  plate 

bending  nultlpller  for  corner  crack  at  a  hole  in  a  plate 

bending  nultlpller  on  etrees  intensity  for  renote  bending 

bending  nultlpller  for  surface  crack  in  a  plate 

half-length  of  cracked  plate 

stress- intensity  factor  (node  I) 

applied  bending  nonent 

curve  fitting  functions  defined  in  text  (i  •  1,  2,  or  3) 
shape  factor  for  elliptical  crack 
radius  of  hole 

renote  bending  stress  on  outer  fiber,  3H/bt* 
renote  unifora  tension  etrees 

thickness  or  one-half  plate  thickness  (see  Fig.  2) 
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X 


function  defined  In  text 


v  FoIiiob'i  ratio  <w  •  0.3) 

4  paraaatrlc  angle  of  ellipse,  deg 


STRESS-INTENSITY  EQUATIONS 

The  stress-intensity  factor,  K,  at  any  point  along  the  crack  front  in  a 
f Inlte-chicknesa  plate,  such  ss  that  shown  In  Figure  1,  was  taken  to  be 


where 

-  (k,  ♦  ♦  H,{f )*],(,.. 


<!•> 


(lb) 


and 


-  Hj  ♦  (H2  -  Bt)  slttP  +  (it) 

The  function  Q  Is  the  shape  factor  for  an  ellipse  and  la  given  by  the  square 
of  the  complete  elliptic  integral  of  the  second  kind  (2].  The  boundary- 
correction  factor,  Fj,  accounts  for  the  Influence  of  various  boundaries  and 
Is  a  function  of  crack  depth,  crack  length,  hole  radius  (where  applicable), 
plate  thickness,  place  width,  sad  Che  paranecrlc  angle  of  the  ellipse.  The 
product  HjFj  ts  Che  corresponding  bending  correction.  The  subscript  J 
denotes  che  crack  configuration:  J  ■  c  is  for  a  corner  crack  in  a  plate, 

)  -  e  la  for  an  eabedded  crack  la  a  plate,  J  ■  a  is  for  a  surface  crack  In 
a  place,  J  •  sh  Is  for  a  surface  crack  at  a  hole  In  a  plate,  and  J  •  ch  is 
for  a  corner  crack  at  a  hole  in  a  plate.  Functions  Np  Mj,  M3,  R^ ,  «2. 

end  p  are  defined  for  each  appropriate  configuration  and  loading.  The 
series  containing  Mt  la  the  boundary-correction  factor  at  the  naxinun  depth 
point.  The  function  f ^  Is  an  angular  function  derived  f roe  che  solution  for 
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an  elliptical  crack  In  an  Infinite  solid.  This  function  accounts  for  aost  of 
the  angular  variation  In  atrees-intensity  factors.  Tha  function  ftf  is  a 
flnlta-vldth  correction  factor.  Function  g  danotas  a  product  of  functions, 
such  as  •••***  chat  are  usad  to  fine-tune  tha  aquations.  Functions 

and  H2  ara  banding  ailtlpliers  obtained  froa  banding  results  at  4  equal 
zero  and  x/2,  respectively.  Figure  3  shoes  the  coordinate  systaa  usad  to 
define  tha  paraaatric  angle,  4,  for  a/c  lass  than  and  a/c  greater  than 
unity. 

Vary  useful  caplrieal  express  Ions  for  Q  hsva  been  developed  by  lave 
(see  Ref.  9).  The  expressions  are 

for  i  <  1  (2a) 

c 

for  ~>  1  (2b) 

c 

The  aaxlsuai  error  In  the  stress-intensity  factor  caused  by  using  these 
approxiaate  equations  for  Q  Is  about  0.13  percent  for  all  values  of  a/c. 
(Rawe's  original  equation  vas  written  In  terns  of  a/2c). 

In  the  following  sections,  the  stress- Intensity  factor  equations  for 
eabedded  elliptical  cracks,  seal-elliptical  surface  cracks,  quarter-elliptical 
corner  cracks,  saxl-elllptlcal  surface  cracks  at  a  hole,  and  quarter- 
elliptical  corner  cracks  at  a  hole  in  finite  plates  (see  Fig.  2)  subjected  to 
either  reaote  tension  or  bending  loads  are  presented.  The  psrtlcular  func¬ 
tions  chosen  were  obtained  froa  curve  fitting  to  f lnlte-eleaent  results  (8,  9, 
18,  and  19]  by  using  polynoalsls  in  terns  of  a/c,  a/t,  sod  angular  func¬ 
tions  of  4.  For  cracks  eaanatiog  froa  holes,  polynoalsl  equations  In  teras 
of  c/r  and  4  vtre  also  used.  Typical  results  will  be  presented  for 


Q  -  1  ♦  l.464(J) 


1.63 


Q  -  l  -r  l-A64(^) 


1.63 
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a/c  ■  0.2,  0.5,  1,  and  2  with  a/t  varying  froa  0  to  l.  Table  1  given  the 
range  of  applicability  of  h  •/*,  a/e,  r/t,  and  (r  +  e)/b  for  the 
propoeed  aquations. 

Babedded  Elliptical  Creek 

The  stress-intensity  factor  equation  for  an  oabedded  elliptical  crack  in 
a  finite  plote.  Figure  2(a),  subjected  to  tenaloa  ana  obtained  by  fitting 
equation  (1)  to  finite-elans at  reaultn  in  Reference  19.  The  reeolte  of 
Irvin  (1)  were  ueed  to  account  for  the  Halting  behavior  aa  a/c  approached 
zero  or  Infinity.  The  equation  la 


K  -  s  fS  '.<!•  ?•  t  ♦> 


(3) 


for  0  <  a/c  <  -,  c/b  <  0.S,  and  -«<♦<*  provided  that  a/t  aatiaflea: 


±  <  1.23  0.6)  for  0  <  i  <  0.2 

C  'C  c 


(*) 


f<* 


for  0.2  <  i  <  • 
c 


The  function  r,  accounts  for  the  Influence  of  crack  shape  (a/c),  crack  else 
(a/t),  finite  width  (c/b),  and  angular  location  (9),  and  was  chosen  as 


'.  •  ["i  ♦  ♦  ")(!*]*  '♦ f.  <5) 

The  tens  in  brackets  gives  the  boundary-correction  factors  at  p  •  s/2  (where 
g  -  f ^  •  l).  The  function  f^  was  taken  froa  the  asset  solution  for  aa 
eabedded  elliptical  crack  in  aa  Infinite  solid  (l)  and  f„  is  a  finite-width 
correction  factor.  The  function  g  is  s  fine-tuning  curve-fitting  function. 
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for  a/c  <  It 


*1 

*4 


1 


o.os 

- 7J7T 

o.n  «■  (£) 


*S 


0.29 


0.23  + 


and 


t 


(£)  (2*«  “  2 

1  ♦  *(£) 

I 


(CM  ♦! 


cos  #  ala 


(6) 

(7) 


(•) 


(9) 


(10) 


(Hoc*  chat  aq.  (9)  la  allghcly  different,  and  la  believed  ta  ha  nor#  accurate, 
chan  that  glean  la  Kaf.  19.)  The  finite-width  correction,  f^,  ftoa 
Reference  9  vea 

for  c/h  <  0.3.  (Race  that  far  the  a^iddad  crack,  t  la  daflaad  aa  one-half 
of  the  full  plate  thlckaaaa.) 


Far  a/c  >  1: 


*1  * 


(12) 
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Th«  functions  Mj,  N),  g,  sad  fw  arc  th*  um  as  aquae  Iona  (7),  (I),  (9), 
and  (11),  respectively. 

Figure  4  shows  soee  typical  boundary-correction  face or a  for  varloua  crack 
shapes  (a/c  *  0.2,  0.5,  1,  and  2)  with  a/c  aqual  to  0,  0.5,  0.75,  mod  l. 

Ttaa  co rr ace Ion  factor,  F#,  la  plot tad  agalnat  tha  paraaacar  angle,  4.  At 
4*0,  eha  point  on  tha  crack  frost  that  la  locatad  at  tha  castor  of  tha 
plata,  tha  lafluasca  of  plata  thlcksaaa  la  aieh  laaa  than  at  4  •  e/2,  tha 
poise  that  la  locatad  cloasae  to  eha  plata  surface.  Tha  raaolta  shows  for 
a/t  ■  0  ara  tha  exact  aolutlooa  for  as  alllptlcal  crack  Is  as  Infinite  solid 
(1].  For  a/t  <  0.8,  the  raaolta  fros  eha  aquation  ara  within  about  3  poreant 
of  eha  flntea-alasant  results.  (Haratn,  “percent*  error  Is  defined  aa  the 
difference  between  the  aquation  and  the  f Inlte-aleneot  raaulta  00 real lead  bp 
eha  eaxleue  value  for  that  particular  case.  This  definition  la  necessary 
because.  In  sons  cases,  tha  streee-’loteesity  factor  r safes  fros  positive  to 
negative  along  the  crack  front.)  For  a/t  >  0.8,  the  accuracy  of  aquation  (3) 
has  not  been  established,  lot  Its  use  In  that  range  appears  to  oa  supported 
by  estlaatee  baaed  on  aa  snbedded  crack  approaching  a  through  crack  (see 
Ref.  19). 

•ending  aquations  were  not  developed  for  tha  esbedded  alllptlcal  crack. 

Sesl-f lllpt teal  Serfaca  Crack 

The  aquations  for  the  strese-lateneity  factors  for  a  sanl-elllptlcal 
surface  crack  in  a  finite  plata,  Figure  2(b),  subjected  to  resote  tension  *nd 
bending  loads  ware  obtained  fros  kef erases  9.  Tha  tension  and  beading 


equations  wrt  previously  fitted  to  flnlte-eleaent  results  fron  Raju  sod 
.tswnen  (8)  for  s/e  values  less  then  or  equal  to  unity.  Equations  fer 
tension  end  heeding  loads  for  s/e  greater  than  unity  were  developed 
herein.  The  results  of  Cross  and  Srawley  ( 20)  for  a  single-edge  crack 
were  used  to  aceouat  for  the  Halting  behavior  aa  a/e  aproaehes  taro. 

The  equation  Is 

*-<*«*  ».v  fi\  f-  f-  ♦)  a*> 

for  0  <  a/c  <2,  e/b  <0.5,  and  0  <  g  <  «,  again,  provided  that  a/t 
satisfies  equation  (4).  The  function  F#  was  chosen  to  be 


r 


e 


(15) 


for  a/c  <  1: 


Hj  •  1.13  -  0.09(2-) 


(16) 


ft,  •  -0.54  f 

^  0.2  ♦  2 

c 


(17) 


*3 


0.5 


1 

0.45  ♦  1 
c 


♦  14(1 
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l  ♦  0.1  ♦  0.35(2)*  (l 


sin  *)' 


(18) 


(19) 


*•6  f j  i*  givet.  bf  equation  (10).  The  fiolte-vidth  correction,  (v,  la 
again  given  by  equation  (11).  Equations  (IS)  through  (19)  were  taken  fron 
Reference  9.  (The  Urge  power  la  eq.  (18)  was  needed  to  fit  the  behavior  as 
a/c  approaches  tero.) 
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The  bending  multiplier ,  Hj,  In  equation  (1)  has  the  fora 

HJ  *  H1  +  (H2  "  Hi.)  8lnP  ♦  (20 

where  Hlt  H2,  and  p  are  defined  for  each  crack  configuration  conaldered. 
For  the  surface  crack  (j  ■  s). 


p  -  0.2  +  -  +  0.6  - 

r  c  t 


a,  •  l  -  0.34  «  -  o.u  |(i) 


»2  -  1  +  °2l(f)  *  °22(t) 


In  this  equation  for  H2, 


G_.  -  -1.22  -  0.12  2 

L  l  C 


0.75  1.5 

C„  -  0.55  -  1.05(i)  *0.47(±) 


Equations  (21)  through  (25)  were  taken  from  Reference  9. 


For  a/c  >  l: 


Ml  “  l/f  U  ♦  0-0«  f) 


*2  "  °'2(f) 


m3  -  -0. 11(f) 


l 

i 


(29) 


g  -  1  + 


jo. I  +  0.35(f)jf  *j(l  - 


sin  $) 


and  and  fw  are  given  by  equations  (13)  and  (11),  respectively. 

The  bending  aultlpller  for  a/c  >  1  Is  also  given  by  equation  (20)  where 


p  -  0.2  +  -  +  0.6  - 

r  at 

“i  ■ 1  ♦  °n ! ♦  °i2  tfi' 
»2  -  1  ♦  °2l  !  *  °22<f)' 
«u  '  -0'04  *  I 


0.75  1.5 

C12  -  0.55  -  1.93(f)  +  1.38(f) 


(30) 

(31) 

(32) 

(33) 

(34) 


Cn  -  -2.11  +  0.77  f 


(35) 


and 


0.75  1.5 

C22  -  0.55  -  0.72(f)  +  0.14(f) 


(36) 


Figures  5  and  6  show  some  typical  boundary-correction  factors  for  various 
surface  crack  shapes  (a/c  “  0.2,  0.5,  1,  and  2)  with  a/t  equal  to  0,  0.5, 
and  1  for  tension  and  bending,  respectively.  For  all  conblnatlons  of  param¬ 
eters  Investigated  and  a/t  <  0.8,  equation  (14)  was  within  ±5  percent  of  the 
finite-element  results  (0.2  <  a/c  <  2)  and  the  single-edge  crack  solution 
(a/c  •  0).  Por  a/t  >  0.8,  the  accuracy  of  equation  (14)  has  net  been 
established.  However,  Its  use  In  that  range  appears  to  be  supported  by 
estimates  based  on  a  surface  crack  approaching  a  through  crack. 
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The  use  of  negative  stress-intensity  factors  in  the  case  of  bending  are 
applicable  only  when  there  Is  sufficient  tension  to  keep  the  crack  surfaces 
open;  that  Is,  the  total  stress-intensity  factor  due  to  combined  tension  and 
bending  must  be  positive. 

Quarter-Elliptical  Corner  Crack 

The  stress-intensity  factor  equations  for  a  quarter-elliptical  corner 
crack  In  a  finite  plate.  Figure  2(c),  subjected  to  tension  and  bending  loads 
were  obtained  by  fitting  equation  (1)  to  the  finite-element  results  In  Refer¬ 
ence  19  for  tension  and  the  results  In  Table  1  for  bending.  The  equation  Is 


*  ■  (st  *  -A)f!  Fc!?  f-  !•  ♦> 


for  0.2  <  a/c  <  2,  a/t  <  1,  and  0  <  $  <  n/2  for  -  <  0.5.  The  function 


Fc  was  chosen  as 


[»L  *  M2 tfl2  * 


g,  f  .  f 
°2  $  w 


For  a/c  <  1: 


-  1.08  -  0.03(|) 


M,  -  -0.44  +  — -0-6.- 

2  0.3  +± 

c 


L 

M3  -  -0.5  +  0.25(|)  +  14.8(1  -  |) 


{l  -  1  +  Jo. 08  +  0.4(5-)J 


(1  -  sin  $)' 
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g-  -  1  + 


jo.08  ♦  0.1S(±)J(l  -  cos  p)3 


and  f  ^  Is  glvsa  by  aquation  (10).  Tho  finite-width  correct Ion,  fv,  «u 
•stlaatsd  herein  by  using  the  single-edge  crack  tension  solution  given  In 
Reference  21  (divided  by  1.12)  and  was 


fw  -  1  -  0.2X  ♦  8.4X2  -  lt.4X3  ♦  27. IX4 


where  X  *  (The  width  correction  froa  Ref.  21  was  divided  by  1.12 

because  tne  front-face  correction  was  already  Included  In  «q.  (38).) 
Equation  (44)  la  restricted  to  c/b  <  0.5. 

As  a/t  approaches  unity,  with  a/c  •  1  and  p  -  0,  the  stress- 
intensity  factor  equation  (eq.  (37))  for  tension  reduces  to 


K  -  St/^  l.llf 


Equation  (45)  Is  within  about  1  percent  of  the  accepted  solution  [21]  for 
c/b  <  0.6. 

The  bending  wltlpller,  Hc,  has  the  fora  given  by  equation  (20). 
Functions  p,  Hlt  «2,  ar.d  C2l  are  given  by  equations  (2l)-(24), 
respectively,  for  a/c  <  l.  The  function  Is 


,  0.75  1.5 

C22  “  0,64  "  l -05(f)  ♦  0.47(1) 


For  a/c  >  1: 


*1  -/|  fl-08  -  0.03  S) 


M2  -  0.375(f) 
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(49) 


M3  -  -0.25(f) 


t  -  1  +  jo. 08  +  0.4(f)J<l  -  sin  *)3 
+  |o.08  +  0.l5(f)J(l  -  cos  *): 


82  “  1 


(50) 


(51) 


and  f  ^  Is  given  by  equation  (13).  The  finite-width  correction  Is,  again, 
given  by  equation  (44). 

The  bendlng-correctlon  factor  Hc  Is,  again,  given  by  equation  (20) 
where  p,  ,  H2,  Gllt  Cl2,  and  G2l  are  given  by  equations  (30)-(35), 

respectively.  The  function  C22  Is  given  by 

0.75  1.5 

G22  «  0.64  -  0.72(f)  +  0.14(f)  (52) 


Figures  7  and  8  show  some  typical  boundary-correction  factors  for  corner 
cracks  In  plates  for  various  crack  shapes  (a/c  -  0.2,  0.5,  l,  and  2)  with 
a/t  varying  from  0  to  1  for  tension  and  bending,  respectively.  At  a/t  **  0, 
the  results  for  tension  and  bending  are  Identical.  As  expected,  for  tension 
the  effects  of  a/t  are  tmich  larger  at  lower  values  of  a/c.  Again,  the  use 
of  negative  stress-intensity  factors  In  this  case  of  bending  are  applicable 
only  when  there  is  sufficient  tension  to  keep  the  crack  surfaces  open  (stress- 
intensity  factor  due  to  combined  tension  and  bending  must  be  positive). 

Seml-Elllptlcal  Surface  Crack  at  Hole 
Two  symmetric  surface  cracks.-  The  stress-intensity  factor  equation  for 
two  symmetric  semi-elliptical  surface  cracks  located  along  the  bore  of  a  hole 
In  a  finite  plate,  Figure  2(d),  subjected  to  tension  was  obtained  by  fitting 
equation  (l)  to  finite-element  results  (19].  The  equation  Is 


15 


ORtGJNAL  PASS  T3 
*)  OF  POOR  QUALfTY 


(53) 


K  "  St  \/*f  F«h^c*  f‘  P  ¥*  ¥• 


Cor  0.2  <  a/c  <  2,  a/t  <  1,  0.5  <  r/c  <2,  (r  +  c)/b  <  0.5,  and 

-it/ 2  <  t  <  it/ 2 .  (Note  Chat  here  t  la  defined  a a  one-half  of  Che  full  place 

chlckness.)  The  funcClon  F  h  was  chosen  as 


-  ["l  +  "2(t)2  *  »2  *3  f. 


For  a/c  <  1: 


■  l 


M2 - 3 

0.11  ♦  (f) 


0.23  +  (i) 


*1  “  1 


-  4  a  1/2 
(f)  (2.6  -  2  f) 

1  ♦  4(?) 


1  +  0.3581  +  1.425X2  -  1.57813  +  2.156X* 
l  +  0.08X2 


l  +  —  cos(0.9$) 


7  M  1 

g3  *  1  +  0.1(1  -  cos  ♦)  fl  -  i) 


(Noce  thac  eq.  (58)  Is  sllghely  different,  and  Is  believed  Co  be  aore  accurate. 
Chan  chac  given  In  Ref.  19.)  The  function  f^  Is  given  by  equation  (10). 


where  n  •  1  le  for  a  single  crack  and  n  •  2  Is  for  two-syaaetrlc  cracks. 
This  equation  was  chosen  to  account  for  the  effects  of  wtdth  on  stress  concen¬ 
tration  at  the  hole  [22]  and  for  crack  eccentricity  [21].  For  a/c  >  1: 


The  functions  M2,  M3,  glt  g2 ,  gj,  and  X  are  given  by  equations  (56) 

through  (61),  and  the  functions  f^  and  ftf  are  given  by  equations  (13)  and 
(62),  respectively. 

Eatlsates  for  a  single-surface  crack.-  The  stress-intensity  factors  for  a 
single-surface  crack  located  along  the  bore  of  a  hole  were  estimated  fro*  the 
present  results  for  two  symaetrlc  surface  cracks  by  using  a  conversion  factor 
developed  by  Shah  [15].  The  relationship  between  one-  and  two-surface  cracks 
was  given  by 


where  K  for  two  cracks  oust  be  evaluated  for  an  Infinite  plate  (fy  »  1)  and 
then  the  finite-width  correction  for  one  crack  Bust  be  applied.  Shah  had 
assumed  that  the  conversion  factor  was  constant  for  all  locations  along  the 
crack  front,  that  Is,  Independent  of  the  parsaetrlc  angle. 

Figure  9  shows  some  typical  bounds ry-corrsctlon  factors  for  a  single 
surface  crack  at  a  hole  for  various  crack  shapes  (a/c  •  0.2,  0.5,  l,  and  2) 
with  a/e  varying  fro*  0  co  l.  These  results  were  in  good  agreeaent  with 
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boundary-correction  factors  est'laati  •:  by  Shah  [IS].  The  agreement  was 
generally  within  about  2  percent  except  where  the  crack  Intersects  the  free 
surface  (2*/n  ■  1).  Here  the  equation  gave  results  that  were  2  to  5  percent 
higher  than  those  est lasted  by  Shah. 

Stress-Intensity  factor  equations  for  bending  were  not  developed  for  a 
surface  crack  located  at  the  center  of  a  hole. 


Quarter-Elliptical  Corner  Crack  at  a  Hole 
Two  syiesetrlc  corner  cracks.-  The  stress- Intensity  factor  equations  for 
two  symmetric  quarter-elliptical  corner  cracks  at  «  hole  In  a  finite  plate. 
Figure  2(e),  subjected  to  remote  tension  and  bending  loads  were  obtained  by 
fitting  to  finite-element  results  In  Reference  18.  The  equation  Is 


+  %V  fJ  b  b  s-  I-  »> 


for  0.2  <  a/c  <  2,  a/t  <  l,  0.5  <  r/t  <2,  (r  +  c)/b  <  0.5,  and 
0  <  o  <  x/2.  The  function  F.h  was  chosen  as 


[MI  *  »2<f)2  *  "j(f)  j  *1 


*2  *3  *»  '♦  f. 


For  a/c  <  1: 


M,  -  1.13  -  0.09(2-) 

I  c 


.  0.89 

M,  •  -0.54  +  - - 

1  0.2  +  - 

c 


M,  •  0.5 - - - -♦  14(1  -  f) 

1  0.65  +  *  c 

c 
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.V 


(4 


*i  *  1  + 


jo.l  +  0.35(|)J(l  -  aln  g)2 


OiJQd^A-  ?;V.?C 

OF  PCCrt  QUALITY 


1  +  0.358X  •f  1.425X2  -  1.578X3  +  2.15614 
1  +  0.13>2 


where 


1  +  J  COS(pg) 

H  •  0.85  for  tension  end  p  “  0.85  -  0.25(a/t)*^*  •  for  bending.  The 
functions  gj  and  g4  are  given  by 

g3  -  fl  +  0.04  i)  [l  ♦  0.1(1  -  cos  g)2]  jo.85  ♦  0.15 (i)  ^ 


g4  -  1  -  0.7 fl  0.2)(1  -  f)  ( 

Functions  f^  and  ftf  are  given  by  equations  (10)  and  (62),  respectively. 

The  bending  multiplier,  Hch,  Is  given  by  equation  (2C)  for  a/c  <  l. 
The  terms  p,  H^,  and  H2  are  given  by 

P  -  0.1  +  1.3  iti.ii*  0.7  -£)  < 

r  t  2  c  vt ' 

"l  -  1  *  cn  F+C.2'f!2+Cuff)3  1 


»2.i.c21i.c22(i)  .Cj,(i)5 


where 


:u  -  -0.43  -  0.74  i  -  0.84 (i) 


Cl2  -  1.25  -  1.19  £  +  4 . 39 fi ) 
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(4 


2 

Cn  -  -1.94  +  4.22  f  -  5.51(f) 

2 

C21  *  -1-5  -  0.04  f  -  1.73(f) 

2 

G22  -  1.71  -  3.17  f  +  6.84(f) 
C23  -  -1.28  +  2.71  f  -  5.22(f)2 


ORIGINAL  PAGE  18 
OF  POOR  QUALITY 


For  a/c  >  1 : 


Mt  -  /f  (l  +  0.04  £) 


M2  -  0.2(f) 


4 

m3  -  -0.11(f) 

gj  -  l  +  jo.l  +  0.35(f)(f)J<i  -  sin  ♦)' 


(SO) 

(81) 

(82) 

(83) 

(84) 

(85) 

(86) 

(87) 


Functions  g2  and  \  art  glvan  by  equations  (71)  and  (72).  Function  g3  la 
given  by 


g3  -  1.13  -  0.09 


8* 


+  0.1<1  -  cos 


♦>2]  0.85  +  0.l5(f)l/4J 


(88) 


and  g&  “  l.  The  functions  f ^  and  fy  are,  again,  given  by  equations  (13) 
and  (62),  respectively. 

Again,  the  bendlng-correctlon  factor,  «ch>  Is  given  by  equation  (20). 
The  function  p  Is  given  by  equation  (30)  for  a/c  >  l.  The  H-functlons  are 
given  by  equatlona  (75)  and  (76)  where 
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T 


I* 


a 


opKssm.  pacu-  m 
OF  POOR  QUALITY 


Gu  •  "2*07  °*06  7  (89) 
Cl2  -  4.35  ♦  0.16  |  (90) 
C13  *  _2-93  '  °-3  f  (91) 
c2l  ’  *3*64  +  °*37  |  (92) 
C22  -  5. 07  -  0.49  |  (93) 
C21  *  "4*32  *  0.53  |  (94) 


Estimates  for  a  a Ingle-corner  crack.-  The  stress-intensity  factors  for  a 
•Ingle-corner  creek  ec  e  hole  were  estimated  from  the  present  results  for  tvo- 
•ymmetrtc  corner  cracks  by  using  the  Shsh-converslon  factor  (eq.  (64)).  Raju 
and  Newman  (18]  have  evaluated  the  use  of  the  conversion  factor  for  some 
corner-crack-at-a-hole  configurations.  The  stress-intensity  factor  obtained 
using  the  conversion  factor  were  In  good  agreement  with  the  results  from  Smith 
and  Kullgren  (16]  tor  a  single-corner  crack  at  a  hole. 

Figures  10  and  11  show  some  typical  boundary-correction  factors  for  a 
single  corner  crack  at  a  hole  for  various  a/c  and  a/t  ratios  for  tension 
and  bending,  respectively.  Again,  the  use  of  negative  stress- Intensity 
factors  in  the  case  of  bending  are  applicable  only  when  there  Is  sufficient 
tension  to  make  the  total  stress- Intensity  factor,  due  to  combined  tension  and 
bending,  positive. 
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CONCLUDING  REMARKS 


Stress-Intensity  factors  fro*  three-dlaensional  f Inlce-eleaent  analyses 
were  used  to  develop  stress- Intensity  fsetor  equations  for  a  vide  variety  of 
crack  configurations  subjected  to  either  reaote  uni for*  tension  or  bending 
loads.  The  follovlng  configurations  eere  Included:  an  embedded  elliptical 
crack,  a  seal-elliptical  surface  crack,  a  quarter-elliptical  corner  crack,  a 
seal-ell lpt leal  surface  crack  along  the  bore  of  a  hole,  and  a  quarter- 
elliptical  corner  crack  at  the  edge  of  a  hole  In  finite  plates.  The  equations 
cover  a  wide  range  of  configuration  paraaeters.  The  ratio  of  crack  depth  to 
place  thickness  (a/c)  ranged  froa  0  to  1,  the  ratio  of  crack  depth  to  crack 
length  (a/c)  ranged  froa  0.2  to  2,  and  the  ratio  of  hole  radius  to  plat* 
chlcknees  (r/c)  ranged  froa  0.5  to  2  (where  applicable).  The  effects  of  plate 
width  (b)  on  acres*- Intensity  variations  along  the  crack  front  were  also 
Included,  but  were  baaed  on  engineering  esc taste*. 

For  all  configurations  for  which  ratios  of  crack  depth  to  plat*  thickness 
do  not  exceed  0.8,  the  equations  are  generally  within  5  percent  of  the  flnlte- 
eleaenc  results,  except  where  the  crack  front  Intersects  a  free  surface.  Here 
the  proposed  equations  give  higher  stress-intensity  factors  than  the  finite- 
element  results,  but  these  higher  values  probably  represent  the  Halting 
behavior  as  the  aesh  Is  refined  near  ehe  free  surface.  For  ratios  greater 
chan  0.8,  no  solutions  are  available  for  direct  comparison;  however,  ch* 
equations  appear  reasonable  on  the  basis  of  engineering  estimates. 

The  stress-intensity  factor  equations  presented  herein  should  be  useful 
for  correlating  and  predicting  fattgue-crack-grovth  races  as  well  as  In 
coaputlng  fracture  toughness  and  rracture  loads  for  these  types  of  crack 
conf lguraclona. 
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Table  1.  Range  of  applicability  for  stress-intensity  factor  equations. 
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<u>a/t  <  1  for  remote  tension  and  a/t  <  0.8  for  remote  bending. 


% 


4 


Table  2.  Boundary-correction  factors,  FCHC,  for  quarter-elliptical 
corner  crack  In  a  place  subjected  to  bending 
(v-  0.3;  FcHc  -  K/(SbvWQ  )). 


a/t 

a/c 

2 

0.2 

0.5 

0.8 

0 

0.522 

0.609 

0.779 

0.25 

0.669 

0.702 

0.808 

0.2 

0.5 

0.801 

0.746 

0.716 

0.75 

0.868 

0.746 

0.577 

1.0 

0.876 

0.750 

0.604 

0 

0.740 

0.799 

0.904 

0.25 

0.724 

0.690 

0.670 

0.4 

0.5 

0.785 

0.632 

0.451 

0.75 

0.826 

0.583 

0.272 

1.0 

0.846 

0.569 

0.262 

0 

1.084 

1.046 

1.027 

0.25 

0.934 

0.770 

0.604 

1.0 

0.5 

0.838 

0.547 

0.237 

0.75 

0.798 

0.417 

0.011 

1.0 

0.839 

0.407 

-0.032 

0 

0.932 

0.811 

0.734 

0.25 

0.851 

0.623 

0.442 

2.0 

0.5 

0.761 

0.413 

0.105 

0.75 

0.700 

0.268 

-0.131 

1.0 

0.677 

0.215 

-0.206 
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Fig.  1 — Corner  cracks  at  Che  edge  of  a  hole  in  a 
finite  plate  subjected  to  remote  tension 
and  bending. 
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Fig.  2 — Embedded-,  surface-,  and  corner-crack  configurations  (all  cracks  have 
elliptical  fronts). 
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20/ 7T 


(C)  0/C  »  1. 


(d>  o/c  -  2. 


Fig.  9— Typical  boundary  correct ion  factors  for  a  single  surface  crack 
at  the  center  of  a  circular  hole  in  a  plate  subjected  to  renote 
tension  (r/t  *  1;  r/b  •  0). 
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o/t  -  0 
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Fig.  n 


(b>  q/c  »  0.5. 


«  v 

2V  * 

(c)  0/C  »  1. 


Q/t  «  0 


(d)  o/c  -  2. 


--Typical  boundary-correction  factor*  for  a  two- syne** trie  corner 
crack  at  the  edge  of  a  circular  hole  In  a  plate  subjected  to 
reaote  bending  (r/t  •  0.5;  r/b  •  0). 


